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Abstract 

We draw motivation from recent experimental studies and present a comprehensive study of 
magnetothermoelectric transport in a graphene monolayer within the linear response regime. We 
employ the modified Kubo formalism developed for thermal transport in a magnetic field. Ther- 
mopower as well as thermal conductivity as a function of the gate voltage of a graphene monolayer 
in the presence of a magnetic field perpendicular to the graphene plane is determined for low mag- 
netic fields (~1 Tesla) as well as high fields (~8 Tesla). We include the effects of screened charged 
impurities on thermal transport. We find good, qualitative as well as quantitative, agreement with 
recent experimental work on the subject. In addition, in order to analyze the effects of modulation, 
which can be induced by various means, on the thermal transport in graphene, we evaluate the 
thermal transport coefficients for a graphene monolayer subjected to a periodic electric modula- 
tion in a magnetic field. The results are presented as a function of the magnetic field and the gate 
voltage. 
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I. INTRODUCTION 



Graphene exhibits remarkable thermal properties. The measured values of thermal con- 
ductivity of graphene reach as high as several thousand of watt per meter Kelvin Q-Q, and 
these are among the highest values of known materials. Heat transport measures the energy 
carried by both electrons and phonons and is fundamental to understanding a material, 
its ground states, excitations and scattering mechanisms. If the dream of carbon-based 
electronics is to be realized, it is essential to study how and how fast heat is dissipated 
across graphene devices. This requires systematic measurements of thermal conductivity 
and thermopower over a broad temperature range (1.5-300 Kelvin) under various external 
conditions. Therefore, recently there has been considerable interest, both experimentaljs)- 
^1 and theoretical 9-16], in the study of thermoelectric and magnetothermoelectric trans- 
port in graphene. This is partly due to the realization that the information provided by 
thermoelectric transport is complementary to electrical transport. And thermoelectric and 
magnetothermoelectric transport studies are extremely useful in providing insight on the 
scattering mechanism involved in transport. Fundamentally related to the electrical con- 
ductivity, the thermal conductivity and thermoelectric coefficients can be determined by 
the band structure and scattering mechanisms. The thermoelectric coefficients involve the 
energy derivatives of the electrical transport counterparts such as the conductivity cxsj. 
Recent measurements of thermoelectric power(TEP) on graphene samples in zero and non- 
zero magnetic fields have shown a linear temperature dependence of TEP which suggest 
that the dominant contribution is that of diffusive thermopower (So)- A comparison be- 
tween the measured TEP and that predicted by the Mott formula shows general agreement, 
particularly at lower temperatures (T < 5070 110|^ However, at higher temperatures devi- 



ation from the Mott relation have been reported f], ij]. In theoretical work, Yan et al. 15] 
have determined the TEP of Dirac fermions in graphene with in the self-consistent Born 
approximation. Also, Hwang et al.[9|, in their calculation of TEP incorporate the energy 
dependence of various transport scattering rates and show that the dominant contribution is 
from the screened charged impurities in graphene's environment. Further, Vaidya et al.[lo| 
used Boltzmann transport theory to calculate Sd in graphene after considering contributions 
of optical phonon and surface roughness scatterings. 

Application of a magnetic field in addition to a thermal gradient has profound effects on 
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the thermal transport in a system and serves as an additional probe. When a magnetic field is 
applied perpendicular to the x—y plane of the sample, the diffusing charge carriers experience 
the Lorentz force. This results in developing a transverse electric field E y in addition to 
the longitudinal field E x . The thermopower is determined from the thermal gradient VT 
and the induced voltage W as S xx = — also known as the Seebeck coefficient) and 
S xy = — -^jKthe Nernst coefficient). They are a measure of the magnitude of the longitudinal 
and transverse voltages generated in response to an applied temperature gradient. They 
are very sensitive in graphene due to its semimetal nature 12]. The quantum magnetic 
oscillations in electrical and thermal transport have been earlier investigated theoretically 
by Gusynin and Sharapov 17| and they obtained analytical results for longitudinal thermal 
conductivity and the Nernst coefficient. However, they assumed a scattering rate that is 
constant in energy, independent of magnetic field and temperature. Hence the self energy 
used is not self consistent. Moreover, they evaluated the longitudinal thermal conductivity 
as a function of the magnetic field at different temperatures but at fixed chemical potential 
and constant impurity broadening. Further, they determined the Nernst coefficient (signal) 
without recourse to the modified Kubo formalism appropriate for thermal transport in a 
magnetic field. They neglected the dependence of T on the chemic al p otential/ carrier 
concentration. Dora and Thalmeier extended the work presented in 17| and studied the 
electric and thermal response of two dimensional Dirac fermions in a quantizing magnetic 
field in the the presence of localized disorder[l8].They evaluated the Seeback coefficient and 
the corresponding thermal conductivity as a function of the chemical potential and the 
magnetic field. They did not determine the Nernst coefficient and the transverse thermal 
conductivity. 

What distinguishes our work on unmodulated graphene from the aforementioned previous 
papers is that we employ the modified Kubo formalism required to study thermal transport 
in a magnetic field. As has been discussed earlier, the usual Kubo formula for thermal 
response functions is invalid in a magnetic field and needs to be modified when calculating 
the transverse (Hall) thermal conductivity and the Nernst coeffecient {"bol. 20]. We use the 



phenomenological transport equations obtained from the modified Kubo formalism |20l. |21 |. 
Further, in the scattering rate and the impurity broadening of the Landau levels the effects 
of the carrier concentration that can be varied by the gate voltage are taken into account. 
In the first stage, we determine the components of magnetoelectrothermal(MET) power and 
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MET conductivity of an unmodulated graphene monolayer in the presence of randomly dis- 
tributed charged impurities. The results are presented as a function of the gate voltage for 
small and large magnetic fields applied perpendicular to the graphene sheet. We determine 
both the Nernst and Seeback coefficients as well as longitudinal and transverse thermal con- 
ductivity. These results are then compared with experimental work. In addition, we have 
also carried out a detailed investigation of the MET transport properties of a graphene mono- 
layer which is modulated by a weak one-dimensional periodic potential in the presence of a 
perpendicular magnetic field. Motivation for this has arisen from recent work, experimental 
and theoretical, that has shown that interaction with a substrate can lead to weak periodic 
modulation of the graphene spectrum. Furthermore, applying patterned gate voltage or 



placing graphene on a pre-patterned substrate can also lead to modulated graphene 
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Placing impurities or adatom deposition can do the same. In a previous work, we have 
computed the electric transport coefficients of electrically modulated graphene 251]. It was 
shown that modulation turns the sharp Landau levels into bands whose width oscillates 
periodically with the magnetic field. This affects the magnetoelectric transport coefficients 
which exhibit commensurabilty (Weiss) oscillations. The origin of these Weiss oscillations 
is the commensurability of the two characteristic length scales of the system: The cyclotron 
diameter at the Fermi energy and the period of the modulation 26]. An interesting feature 
of electronic conduction in the modulated system is the opening of the diffusive (band) 
transport channel in addition to hopping (collisional) transport. Both these contributions 
to MET transport are taken into account in this work. 

In the following section, we present the general formulation of the magnetoelectrothermal 
transport problem and perform the calculation of the thermopower and the thermal conduc- 
tivity of unmodulated graphene as well as graphene subjected to one-dimensional (ID) weak 
periodic modulation. The results for the transport coefficients as a function of gate voltage 
(V^) for unmodulated graphene are discussed in Section III, where we also make a com- 
parison with experimental results. Following this in Section IV, the results for modulated 
graphene as a function of the gate voltage and the external magnetic field are presented. 
The present paper ends with a summary and conclusions. 
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II. THERMAL MAGNETOTRANSPORT COEFFICIENTS 



As mentioned in the introduction, corrections to the usual Kubo formula for transport 
have to be made when studying thermal transport in a magnetic field. This was carried 
out by Luttinger, Smerka, Streda and Oji 
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21] . We employ the modified Kubo formalism 



to determine the thermal transport coefficients from the electrical J e and thermal (energy) 
current densities Jq 
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Here r\ — r\ — e<f> with n the chemical potential, <fi the scalar potential, e the electronic charge 
and T the temperature of the system. The electrical and thermal transport coefficients: the 



electrical conductivity a, thermopower S anc 
from the above expressions, following 
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the thermal conductivity k can be obtained 
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4^0 = 0,1,2) are, in general, tensors where a,u = x, y. These phenomenological transport 
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27J 6${B) = £)%{-B). a^{E) is the zero- 
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[exp(-^r + 1)] 1 is the Fermi Dirac distribution 



coefficients satisfy the Onsager relation 
temperature conductivity and f(E) = 
function with r\ the chemical potential. The quantity p^ u = is the resistivity tensor 

whose components are p xx = cr yy /A, p yy = a xx /A, p xy = -p yx = a yx /A with A = a xx o yy - 

In order to calculate the thermal transport coefficients for graphene, we consider a 
graphene monolayer in the xy — plane subjected to a magnetic field B along the z-direction. 
In the Landau gauge, the unperturbed single particle Dirac-like Hamiltonian may be written 

as 

H = v F a. (-ihV + eA) . (7) 



m nMy = 4= ( ~Zu„ Z w/i ) . (8) 



Here, a = {<J x ,a y } are the Pauli matrices and vp = 10 6 m/s characterizes the electron 
velocity with A = (0, Bx, 0) the vector potential. The normalized eigenfunctions of the 
Hamiltonian given in Eq. (f7|) are 

Jk y y /_,0 n „ 1 [( X + Xo )/iy 
y/2L~l\ (Pn[(x + X )/1\ 

where <f> n (x) and n _i(x) are the harmonic oscillator wavefunctions centred at x Q = l 2 k y . 
n is the Landau level index, I = the magnetic length and L y the length of 2D 

graphene system in the ^/-direction. The corresponding eigenvalue is E n = fhj g y/n where 
u g — v f ^/2eB/h = vpy/2/l is the cyclotron frequency of the Dirac electrons in graphene. 

In order to investigate the effects of modulation, we express the Hamiltonian in the 
presence of modulation as H = H Q + U(x). Here, U(x) is the one-dimensional periodic 
modulation potential along the x-axis. It is given by U(x) = V e cosKx such that K = — , 
a is the period of modulation and V e is the constant modulation amplitude. To account for 
weak modulation, we take V e to be an order of magnitude smaller than the Fermi energy 
Ep = vpHkp,where kp = y/2irn e is the magnitude of Fermi wave vector with n e the density 
of electrons. This allows us to apply standard first order perturbation theory to determine 
the energy eigenvalues in the presence of modulation. Thus, energy eigenvalues for weak 
modulation (V e < Ep), are E ntky = E n + F n)B cosKx. Here, F^ B = y exp(-| )[L n (u) + 
L n -i(u)), u = and, L n {u) and L n -%(u) are Laguerre polynomials. 

In the presence of a periodic modulation, there are two contributions to magnetocon- 
ductivity: the collisional (hopping) contribution and the diffusive (band) contribution. The 
former is the localized state contribution which carries the effects of Shubnikov de Hass 
(SdH) oscillations that are modified by periodic modulation. The diffusive contribution is 
the extended state contribution and arises due to finite drift velocity acquired by the charge 
carriers in the presence of modulation. In the linear response regime, the conductivity tensor 
is a sum of a diagonal and a non diagonal part : cr^ v (uj) = cr*(u) + a^(cu), fi, v — x, y. In 
general, a^iuS) = a^J^co) + a^(cj), accounting for both diffusive and collisional contribu- 
tion whereas cr^(o;) is the Hall contribution. Here, a xx = ex™ 1 and a yy = ^ x ° x l + (y y h J^. Similar 
to the conductivity tensors, the diagonal components of the thermal transport coefficients 
are determined by the following expressions: 

f(a) _ r-(a)co\ _ p(a) col /q\ 
'-'xx '-'xx yy y J 
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The finite temperature conductivity components a^ v have been evaluated in [25( for scatter- 
ing by random screened charged impurities of density Nj with impurity broadening T. The 
screened potential (in Fourier space) is U Q = 2ne 2 / e^Jq 2 + k 2 , which is valid for small wave 
vectors, q <C k s , k s being the inverse screening length and e the dielectric constant. There- 
fore, from Eq.(|6]), we obtain the zero-temperature phenomenological transport coefficients 
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where r is the scattering time. Here, we have taken the scattering time to be independent 
of Landau-level index n. And the components of thermopower are given by the following 
equations: 
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The components of the thermal conductivity are given by 
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III. MAGNETOTHERMOPOWER AND MAGNETOTHERMAL CONDUCTIV- 
ITY OF UNMODULATED GRAPHENE 



From the electrical conductivity 07^, calculated in our previous work [25|, we determine 



the phenomenological transport coefficients C y a x employing EqsJTTj [T2l and [T31 Employing 



ents jL^yz 

these, the components of thermopower and thermal conductivity are numerically evaluated 
using EqslMl through I2T1 The results for the magnetoelectrothermal transport properties of 
an unmodulated graphene monolayer as a function of the gate voltage are presented in this 
section. The number density n e is related to the gate voltage V g through the relationship n e = 
e eV g /te, where e G and e = 3.9 are the permittivities for free space and the dielectric constant 
for graphene on a SiC>2 substrate, respectively. The electron charge is e and t(~ 300nm) 



is the thickness of the sample [30j. The components of thermopower (S^) and thermal 
conductivity (n^u), as the system moves away from the charge neutral point on the electron 
side on changing the gate voltage, are shown in Fig. (1) at a magnetic field of one Tesla. 



The lattice temperature of 10K and mobility of jj, = 20m /Vs^^ is chosen. The scattering 



time is related to the mobility as r = in a graphene monolayer 32] . Since impurity 



broadening T can be ex pre ssed in terms of the self energy E (E) as T = T(E) = 2 Im [E (E)] 
and also T(E) = h/r 



331 ] . We use the expression for Im [E (E)] derived in |25j to find 



T = ^H(huj g ) 2 1 (AtctEf). The electron number density is n e = 7.19V g x 10 14 m 2 and Fermi 



energy is Ep = hvF\/2irn P . = AA.Sy/VgmeV . And the impurity density is related to V through 



Nj = iil 2 T 2 /Ug 34j. The scattering time of r = 4.431/Xa/T^ x 10 14 s, impurity broadening 



T = 5.93 4 -\/ B / { iiVg)meV and impurity density Nj = x 10 14 m -2 were employed in this 
work 3lM37j|. Moreover, the same study is carried out at a higher magnetic field of 8.8T for 
graphene with mobilities of fi = lm 2 /Vs and /i = 20m 2 /Vs respectively and the results are 
shown in Fig. (2). Since S xx and S yy are identical so only S xx is depicted in these figures. 
The longitudinal coefficient of thermopower (S xx ) is equivalent to the Seebeck coefficient 
and our results provide a qualitative as well as quantitative understanding of the overall 
behavior of the observed S xx (V g ). S xx can have either sign and it is negative in our case 
since the charge carriers are electrons in this range of V g . The transverse component of 
thermopower (S yx ) is also known as the Nernst signal and it arises due to the presence of 
the perpendicular magnetic field as the Lorentz force bends the trajectories of the thermally 
diffusing carriers. It can be seen from Fig. (la) and Fig. (2a) that S xx follows 1/ \fV~g (with 
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Vg oc n e ). Similar behavior of S xx is observed in experiments 0-7]. Notice that we have 
presented results for diffusive thermopower and we have ignored thephonon contribution 
to thermopower due to weak electron phonon coupling in graphene p, |9j. S xx and S yx ( 
S yx — —S xy ) both show Shubnikov-de Haas (SdH) type oscillations in the Landau quantizing 
magnetic field. At the lower magnetic field of 1 Tesla (Fig. (1)), the oscillations are more 
closely spaced since the separation between the Landau levels, which is proportional to the 
magnetic field strength, is smaller compared to the results for the higher magnetic field 
of 8.8 Tesla, (Fig. (2)). Moreover, we observe in Figs, (la) and (2a) that both S xx and 
S yx approach zero at those values of V g where there are boundaries of Landau Levels and 
no carriers are available to participate in transport. The peaks of S xx are observed at the 
centre of Landau levels. With the increase in V g and hence an increase in n e , higher Landau 
levels are occupied. The oscillations in S xx and S xy are damped as we increase V g . The 
reason for this is that higher V g corresponds to higher values of the Fermi energy and if the 
Fermi energy is much larger that the Landau level separation, Landau quantization effects 
are lost. At B — 8.8T, (Fig. (2a)), the oscillations in S xx and S xy show that the width 
of the peaks broaden compared to those for smaller magnetic field of B = IT. At lower 
magnetic field, the separation between the Landau levels is smaller compared to higher 
fields with the result that the peaks of S xx are more closely spaced. Furthermore, the overall 
magnitude of S xx and S xy increases with increasing magnetic field strength (See Fig. la 
and Fig. 2a). In these figures, we also present thermal conductivity as a function of the gate 
voltage. The longitudinal thermal conductivity k xx shows oscillating behavior which damps 
out as V g increases, where Landau quantization effects become less significant. However, 
the transverse component of thermal conductivity n yx rises monotonically with V g as shown 
in Fig. (lb) and Fig. (2b). At the higher magnetic field, quantum Hall steps have begun 
to appear. The behavior of longitudinal and transverse thermal conductivity follows that 
of the corresponding components of electrical conductivity. At higher magnetic fields, the 
splitting of the peaks in the longitudinal thermal condutivity k xx is seen in Fig. (2b) which 



was also observed in 



18j where it is shown that the splitting occurs in such a way that they 



produce antiphase oscillations with respect to the electric one and lead to the violation of 
the Wiedemann-Franz law. For the un-modulated case, o yy — o xx , jC> xx — ^yy c^nd using 
Eq.t flij) through Eq.( l2"T|) we find that S yx = —S xy , k xx = K yy and n xy = —K yx . Therefore, 
only k xx and K yx are shown in the figures. We find that the results for magnetothermal 
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power obtained in our work at B = 8.8T with T = 10K are in good agreement, both 
qualitative and quantitative, with the experimental results obtained in 0,13], see Fig (3) 
of js]. These results do indicate that scattering from screened charged impurities is the 
dominant scattering mechanism required to explain the experimental results. We must add 
that our quantitative results for S yx depend strongly on the mobility of the graphene system. 



IV. MAGNETOTHERMOPOWER AND MAGNETOTHERMAL CONDUCTIV- 
ITY OF PERIODICALLY MODULATED GRAPHENE 

Now we consider the effects of modulation. The ID modulation broadens the sharp Lan- 
dau levels into bands and gives rise to an additional diffusive (or band) contribution to 
transport. This additional contribution is absent without modulation. We now focus on 
the modulation induced changes in the thermal magnetotransport coefficients of graphene. 
Therefore, in the first part we present the thermopower and thermal conductivity of mod- 
ulated graphene with mobility of 20m 2 /Vs as a function of the gate voltage. These are 
shown in Fig. (3) and Fig. (5) respectively. The results are for a constant external magnetic 
field of B = IT applied perpendicular to the graphene sheet, with electric modulation of 
strength V e = 3meV applied in the x-direction at a temperature of T = 10K. In this case 
T = -^LmeV and Hu g = 36.3meV, such that f < F c < fvjj„ to satisfy the requirements of 



weak modulation. The period of modulation is a = 382nm. The results for S xx and S yy are 
identical, so only S xx is shown in these figures. The amplitude of oscillations in S xx (AS XX ) 
is greater than that of S xy (AS xy ) which damps out with increasing gate voltage(V^). Both 
S xx and S xy show SdH-type oscillations and it verifies that the system is Landau quantized. 
The modulations effects are apparent in S xx and S xy which shows modulation of SdH-type 
oscillations and AS XX 3> AS xy , Fig. (3). K yx is greater than k xx and n yy as shown in Fig. 
(5). These modulation induced effects on thermal transport coefficients can be highlighted 
by calculating the difference between the modulated case and the un-modulated case. The 
contribution of modulation to thermopower AS^ y {y e ) = AS^ u {y e ) — AS^(O) and thermal 
conductivity An^iVg) = AK^ u (V e ) — Ak^ u (0) are shown in Fig. (4) and Fig. (6) respectively. 
These figures clearly show the modulation of SdH oscillations in both the thermopower and 
the thermal conductivity. For an un-modulated however for modulated 

graphene k xx ^ K yy and this expected behaviour is seen in Fig. (6) where An xx ^ An yy . 
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The ID modulation gives a positive contribution to Afi ra while Ak xx and An yx oscillate 
around zero. An yy 3> Ak xx , which is a consequence of the fact that Ak xx has only colli- 
sional contribution, whereas Ak w in addition to the collisional part, has large contribution 
from band conduction. 

We also show the the results when the magnetic field is varied and the electron density 
is fixed at n e = 3.16 x 10 15 m~ 2 which corresponds to a gate voltage of V g = 4.39V. The 



Fermi energy of the system is Ep = hvFy/2irn e w 92.3meV. We have taken the mobility of 
20m 2 /Vs 3l|and hence scattering time is taken to be r = 1.86 x 10~ 12 s. Impurity broadening 
T = O.G33\/BmeV and impurity density Nj = 1.23 x 10 13 m~ 2 were employed in this part 
of the work. The strength of the electrical modulation is taken to be K = 2meV with 
period a = 382nm and temperature T = 10K. The difference between the modulated case 
and the unmodulated case highlights the modulation induced effects in these thermoelectric 
quantities. The thermopower and the change in thermopower due to modulation AS fJiV {B) 
are shown in Fig. (7) as a function of the magnetic field in the units of — fce/e. When B 
is less than 0.2T Weiss oscillations are observed whereas SdH type oscillations dominate at 
higher magnetic fields. It is also seen that these oscillations in S xx are 90° out of phase 
with those in S xy . The amplitude of the oscillations in ^> AS xy and they are 90° 

out of phase. Again for B greater than 0.2T the oscillations appear as envelopes of SdH 
oscillations. The different components of the thermal conductivity tensor and the correction 
to it due to ID modulation are shown in Fig. (8). The magnetic field dependence of the 
thermal conductivity tensor is similar to that of the electrical conductivity tensor obtained 
in 25)]. In Fig. (8) we see that An yx ^> Ak vv ^> Ak xx , such that An yx and Ak xx are 180° 
out of phase from each other. 

To conclude, in this work we have studied magnetothermoelectric transport in graphene 
in the linear response regime using the modified Kubo formalism appropriate for thermal 
transport in a magnetic field. Results are presented for both unmodulated graphene as well 
as graphene that is weakly modulated by an electric modulation. We take into account 
scattering from screened charged impurities and our results indicate that these provide the 
most dominant scattering mechanism at low temperatures. The thermopower, the Seebeck 
coefficient and the Nernst coefficient are determined as a function of the gate voltage. Fur- 
thermore, we also determine the magnetothermal conductivity tensor, both the longitudinal 
and the transverse (Hall) components. For unmodulated graphene we were able to make a 
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comparison of the thermopower with experimental results and find that they are in good 
agreement, both qualitative as well as quantitative, with experimental results. In the case of 
modulated graphene, we focus on the modulation induced effects that appear as commensu- 
rabilty (Weiss-type) oscillations in the magnetothermoelectric coefficients. The results are 
presented as both functions of the gate voltage and the magnetic field. 
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